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A mechanism for orientation of bilayer classical Wigner crystals in a piezoelectric medium is 
considered. For the GaAs system the piezoelectric correction to the electrostatic interaction between 
electrons is calculated. It is shown that taking into account the correction due to the piezoelectric 
effect leads to a dependence of the total energy of the electron crystal on its orientation with respect 
to the crystallographic axes of the surrounding matrix. A generalization of Ewald's method is 
obtained for calculating the anisotropic interaction between electrons in a Wigner crystal. The 
method is used to calculate the energy of bilayer Wigner crystals in electron layers parallel to the 
crystallographic planes (001), (0—11), and (111) as a function of their orientation and the distance 
O^l . between layers, and the energetically most favorable orientation for all types of electron lattices in a 

bilayer system is found. It is shown that phase transitions between structures with different lattice 
symmetry in a Wigner crystal can be accompanied by a change of its orientation. 
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I. INTRODUCTION 



It is known that a system of electrons in the presence of a neutralizing positive background at sufficiently low 
temperatures and densities goes into a Wigner crystal state. In particular, this situation takes place in two-dimensional 
electron layers in AlGaAs-GaAs heterojunctions. The formation of a Wigner crystal phase occurs under the condition 
that the average distance between electrons is much greater than the effective Bohr radius (in the absence of an external 
magnetic field) or the cyclotron radius (in high magnetic fields). The latter situation corresponds to a filling factor 
i/<1. In quantum Hall systems the formation of modulated electron structures can occur in other regimes as well. 
For example, skyrmion lattices can arise in the quantum Hall ferromagnet regime {i> « 1, 1/3). jjj In low magnetic 
fields [y Ki N+l/2, where N is an integer and N > 4) the formation of stripe structures can occur at the top, partially 
filled Landau level. The formation of such structures was predicted theoretically || and confirmed experimentally ||] 
by observation of a strong anisotropy of the conductivity. 

An interesting question is that of the orientation of the nonuniform electron structures relative to the crystallographic 
axes of the surrounding matrix. In the two-dimensional electron layers realized in AlGaAs heterostructures, an 
£^ ■ important influence on the orientation of the electron lattice can be exerted by the piezoelectric interaction between 
the electron and elastic subsystems. This possibility was first pointed out in relation to Wigner crystals in Refs. 
and H . A piezoelectric mechanism for the orientation of the stripe structure was considered in Ref . m . In particular, 
in that paper an effect was detected wherein a reorientation of the stripes in bilayer systems arises in the case when the 
period of the structure is greater than the distance between layers. In Ref. || the electron subsystem was described 
using a charge-density-wave model. The purpose of the present paper is to examine the piezoelectric mechanism 
of orientation of nonuniform electronic structures in bilayer systems in another limiting case, corresponding not to 
a charge density wave but to a classical Wigner crystal. As in Ref. ||, we use a model which takes into account 
the anisotropy of the elastic constants of the crystalline matrix. We develop an approach whereby the energy of a 
> system with an anisotropic interaction between electrons can be calculated exactly; this approach is a generalization 
of Ewald's method for calculating Coulomb sums. As a particular case we obtain results pertaining to a monolayer 
system. This topic was discussed previously in Ref. Q , where consideration was limited to the model situation of an 
isotropic crystal. Such a model does not answer the question of the specific orientation that will be realized in the 
GaAs system, the elastic properties of which are described by three elastic constants rather than two. Furthermore, 
the method of rapidly convergent lattice sums was not used in Ref. For the case of bilayer systems, as far as we 
know, the question of a piezoelectric mechanism of orientation of a Wigner crystal has not been considered before. 

The lattice symmetry of a classical Wigner crystal is determined by the minimum of its Coulomb energy. In a 
monolayer system the minimum is achieved for a hexagonal lattice H. In a bilayer system with the same electron 
density in the two layers, five types of electron lattice can form. The structure, dynamical properties, and melting 
criterion of such systems have been studied in detail in Refs. [p|-p0[. Quantum bilayer Wigner crystals in an external 



The possibility of formation of bilayer Wigner crystals in superfiuid 



magnetic field were treated in Refs. |0 and | 1^ 



helium films was examined in Refs. |13[ and pT|. 

In a classical bilayer Wigner crystal the transition between different crystalline phases is regulated by the parameters 
7] = dy/n, where d is the distance between layers and n is the electron density in the layer. The case rj = corresponds 
to a monolayer system with a doubled density, and the case r\ = oo to a system of two noninteracting layers. In both 



1 



cases the minimum energy corresponds to hexagonal lattices, with a period differing by a factor of Therefore, for 
finite r\ there must be transition phases: rectangular, square, and rhombic. 

The piezoelectric interaction, generally speaking, can lead to a shift of the boundaries between phases. In GaAs 
the piezoelectric interaction is rather weak, and so the indicated effect will be small. Nevertheless, since the Coulomb 
interaction in a system with cubic lattice symmetry is isotropic, the piezoelectric correction to the interaction between 
electrons can be important for determining the orientation of the electron crystal. 



II. ENERGY OF THE PIEZOELECTRIC INTERACTION BETWEEN ELECTRONS IN A WIGNER 

CRYSTAL 

Consider an infinite piezoelectric medium. The electrostatic potential ip of an electron placed at the origin of the 
coordinate system is given by the following system of equations: 

divD = 47re<5(r), 

^=0, (1) 

where 



dtp 

Di = -e lk - ^Pi,ki u ki - (2) 



" dx k 

is the electric displacement vector, and 

dtp 

CTik = XiklmUlm - 0l,ik-^ — - (3) 

dxi 

is the stress tensor. Here Ei k is the dielectric tensor, \i k im is the tensor of elastic constants, f3i^k is the tensor of 
piezoelectric moduli, and Uik is the strain tensor. After a transformation to the Fourier components, the system of 
equations (1) reduces to an algebraic system, and one can easily write the solution for the electrostatic potential. Let 
us write it out explicitly for a cubic system, the properties of which are determined by three clastic constants en, C12, 
and C44, one piezoelectric modulus ei 4 , and the dielectric constant e: 

^ = S ~ ^ X 9 v5(t)|(q)iq) + ° iX% W 

where \ — ef 4 /ecn is a small parameter in which the expansion is done, Si(q) is the velocity of sound with polarization 
i in the direction q, and p is the density of the medium. The function P is a homogeneous 8th-degree polynomial of 
the form 

P(q x ,q y , q z ) = q 2 {a iq lq 2 y ql + a 2 ^ ^Sfc) + «3 ^ ^fe, ( 5 ) 

z^fc i^k 

where l,k — x, y, z, and 

ai = cn{2c\ 2 - 2cnci2 + c 44 - 2cnc 44 ), 

2 

a 2 — c n c 44 , 

a.3 = 2 C "( c n + c i2)( c n - c i2 - 2c 44 ). (6) 

The x,y,z axes are directed along the fourfold axes of the crystalline matrix. As we see from formula (4), the 
electrostatic potential contains a correction <5(/? q , the presence of which is due to the piezoelectric interaction. In an 
isotropic crystal, in which the sound velocity is independent of the direction and the coefficient 03 in (5) is zero, the 
correction linear in \ can be represented as an expansion in a finite number of spherical harmonics: 

^q = -^X^ A nmY nm (Q^->Pci), (7) 

n m— — n 

where n is even and n < 6. In Eq. (7) ipq and O q are the polar and azimuthal angles, respectively. The coefficients A nm 
are expressed in terms of the longitudinal and transverse sound velocities. Expansion (7) in the case of an anisotropic 
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crystal will also contain higher harmonics allowed by the symmetry of the system. The coefficients A nm for the 
anisotropic case can be found numerically. A calculation of these coefficients for the GaAs system (en = 12.3 x 10 11 
dyn/cm 2 , en — 5.7 x 10 11 dyn/cm 2 , C44 = 6.0 x 10 11 dyn/cm 2 ) shows that the main contribution to the sum in (7) is 
given by the same harmonics (n < 6) as in the isotropic system (the coefficients of the higher harmonics are at least 
an order of magnitude smaller). It should be noted that, because the relations between the expansion coefficients are 
determined by three rather than two elastic constants, the problem cannot be reduced to isotropic even by limiting 
consideration to the lowest harmonics. 

Using expansion (7), one can easily find the correction to the electrostatic interaction between electrons in a 
piezoelectric medium. After doing the inverse Fourier transformation we find 

2 

5F(r) = -^G(e r ,Vr), (8) 
er 

where 

G(G r , V> r ) = 4tt ^ A nm (- 1) » 2w • m2 Y nm (Or, ■ (9) 

The form of the function G calculated for the GaAs crystal is shown in Fig. [I] (the coordinate axes are chosen along 
the fourfold axes). As we see from Eq. (8), the interaction between electrons contains a contribution corresponding 
to attraction, the strength of which depends on the direction of r and decays as 1/r. Since this decay is as slow as 
that of the Coulomb interaction, the convergence of the lattice sums will be slow. Therefore, for correct calculation 
of the piezoelectric correction to the energy of a Wigner crystal it is better use a method of rapidly convergent sums, 
modified to take into account the anisotropy of the interaction. 




FIG. 1. The anisotropy of piezoelectric interaction between electrons in a GaAs host matrix. 
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Consider a bilayer electron system placed in an infinite piezoelectric medium and oriented normal to a certain 
crystalline direction. In this case it is convenient to make a change of variable in Eq. (9) to the new angles 6 
and ip, measured from the normal to the surface of the electron layer and from an axis lying in the plane of the 
layer, respectively (the overall structure of expression (9) is preserved — only the values of the expansion coefficients 
are changed). If the structure of the electron crystal is assumed fixed, then its orientation is determined by the 
contribution of the terms in (9) which depend on the angle ip r in the chosen reference frame. The latter correspond 
to harmonics with m/0. For calculating the lattice sums we rewrite the dependence of the i/> r -dependent part of the 
interaction between electrons in the form 

2 

V an (r,G r ,ip r ) = -^Yl 2 5 ;m cos z e r sin M e r e""^, (10) 

l>0 \m\>0 

Using the explicit form of the spherical harmonics, we can express the coefficients Bi m in terms of the coefficients 
A nm . Since the index n in Eqs. (7) and (9) takes on only even values, the coefficients Bi m are nonzero only for I and 
m having the same parity. If we keep a finite number of spherical harmonics in expansion (9), then the sum in (10) 
will also contain a finite number of terms. We note that to preserve the point symmetry with respect to the angle 
ip when a finite number of harmonics is taken into account, the transformation to the indicated coordinate system 
should be done in Eq. (4) and then the values of the expansion coefficients in (7) should be found numerically. 

Taking Eq. (10) into account, we write the anisotropic contribution to the energy of the electron crystal in the form 



E an = E? n + E%? + E™, (11) 



where 



..2 



-s-VEi-EpiV" (12) 

|m|>0 R^R' 



corresponds to the contribution of the interaction within the layers, 

e: 



■out _ix v v o v d'|R-R'-c|H 

«" - £ 2^ 2^ ° l ™ 2^ NR-R'-cP + ^KM+'+i)/ 2 1 ' 

i>0 |m|>0 r,r' 11 1 J 

describes the contribution of the interlayer interaction, and E^ gives the correction to the interaction with the 
positive neutralizing background. The vectors R and R' in Eqs. (12) and (13) are lattice vectors, and the vector 
c specifies the displacement of the upper sublattice relative to the lower. The lattice sums in (12) and (13) can be 
reduced to a rapidly convergent form with the use of a modified version of Ewald's method (see Appendix). Employing 
this method gives 

Ne 2 \ 

E an = ±yft{S in + S out ), (14) 

£ 

where N is the total number of particles in the layer, and 

Sin - E B °™{E e im ^$(m,TmR 2 )+iW £ e im ^<S>(m, (15) 

m R/0 G^O 



if I R -I- r 1 1 m I 

S out - J2 B ^{H r , 9 ^TT ^ R+C * (I + H, ™[|R + c| 2 + rf 2 ]) 

lm R [|R+c| 2 +rf 2 ] 2 

e- iGc+im ^ G V{l,m,-^,Trnd 2 )} (16) 
where G are reciprocal lattice vectors. In Eqs. (15) and (16) we have introduced the functions 

^m,x) = x P {lj ~^~_\ X) (17) 



4 



i r i N{hm) 

*{l,m,x,y) = +|m|+ - C N(i,m)+s(xy) 



p/ lj-|TO|i-X S Z / 

1 V 2 J s=0 

[ e - 2 ^( s ,V^- Vy) + (-l)^^eV^F( S , + Vv)] (18) 



In Eq. (18) the are binomial coefficients, N(l,m) — max[(|m| — 0/2, (I — \ m \ — 2)/2], and 

F(s, z) = r(s + 1/2) - sgn(z) 7 (s + 1/2, z 2 ). (19) 

In Eqs. (17)— (19) T(x) is the gamma function, and T(k,x) and "/(k,x) are incomplete gamma functions. We note 
that for I, m equal to zero (these terms have not been taken into account, since they give a direction-independent 
correction to the interaction) the sums (15) and (16) reduce to the known expressions for the isotropic case. 



III. ORIENTATION OF A BILAYER WIGNER CRYSTAL IN A GaAs MATRIX 



Let us use the results of the previous Section to determine the orientation of bilayer Wigner crystals lying in the 
planes (001), (0—11), and (111) in the GaAs matrix. In describing the piezoelectric interaction with allowance for 
anisotropy of the elastic constants, we keep in expansion (7) only the harmonics with n < 18 and \m\ < 12. We note 
that in the cases considered below, harmonics with n > 6 only influence the orientation of hexagonal structures in 
layers parallel to the (001) plane and that of square structures in layers parallel to the (111) plane. In these cases the 
symmetry of the system leads to the vanishing of the contribution of the lower harmonics to the energy of the Wigner 
crystal. 

The structure of a Wigner crystal in a bilayer system is specified by the primitive lattice vectors Ri and R2 and 
the vector c of the relative displacement of the sublattices in adjacent layers. The values of these vectors for the five 
types of lattice considered are listed in Table |. The values of r\ for which the change of lattice symmetry occurs 



were obtained in Ref. 
a rhombic lattice), we 



Since we will have need of the functions 8(rj) (for a rectangular lattice) and ot(rj) (for 
rave repeated the calculations of Ref. ||. According to the results of our calculations, the 
transition between the rectangular and square phases occurs at 77 ~ 0.263, that between the square and rhombic 
phases at 77 w 0.621, and that between the rhombic and double hexagonal at 77 w 0.732. The first two are second-order 
transitions, and the last is first-order. These results reproduce those of Ref. ||. (We will not analyze the transition 
between the simple hexagonal and rectangular phases, which, according to Ref. occurs at r) — 0.006, since at this 
transition the orientation of the electron crystal changes only slightly.) 



TABLE I. Structure of double layer Wi 
primitive vectors of the reciprocal lattice; 
layer electron density. 


gner crystals. Ri, R2 are the primitive vectors of the direct lattice; Gi, 
c, the vector of relative displacement of sublattices in adjacent layers; n, 


G 2 , the 
the one 


Lattice type 


Ri 


R 2 




G 2 


c 


Variable parameter 


n 


simple hexagonal 


(a,0) 


(0, v^a) 






|(1,V3) 




1 

a 2 VZ 


rectangular 


(a, 0) 


(0,aS) 


e.o) 


(o.fl) 


|(1,«5) 


1< S < V3 


1 

a - 6 


square 


(a, 0) 


(0,a) 


e,o) 


(0,2?) 


1(1,1) 




1 


rhombic 


(a,0) 


a(cos a, sin a) 


^(1,-cota) 


(0, 2 * ) 

> ' a sin a * 


| (1 + cos a, sin a) 


I <«<f 


1 

sin ql 


double hexagonal 


(o,0) 


f a a\/3 \ 
( 2 ' 2 > 


— (1 - — ) 




2(1 4=1 




2 

a 2 V3 
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The anisotropic contribution to the energy of a bilayer electron crystal as a function of its orientation in the (001), 
(0—11), and (111) planes is shown in Figs, [?], and Q for various values of the parameter 77 (the energy is given as 
the energy per electron in units of e 2 \y/n/2e). 

For a two-dimensional Wigner crystal lying in the (001) plane, the anisotropic contribution to the energy of the 
simple hexagonal structure does not exceed 2 x 10~ 2 (in the indicated units). A minimum of the energy is reached 
at an angle between Ri and the [100] axis which is a multiple of 30°. The absolute minimum of the energy of the 
rectangular phase occurs when one of the primitive lattice vectors is directed along the [100] axis. At small 77 there 
are also local minima corresponding to angles (3 w ±30° between Ri and one of the fourfold axes. With increasing r\ 
the local minima vanish, and the anisotropy energy increases substantially. A minimum of the energy of the square 
phase corresponds to an orientation of the elementary lattice vectors along the fourfold axes. The vector c here is 
oriented along one of the twofold axes. This orientation of the vector c is preserved in the rhombic phase as well 
(the vectors Ri and R2 rotate smoothly on changes in 77). The transition to the hexagonal phase as the parameter 
77 is increased further is accompanied by a jumplikc change in the orientation of the vectors Ri and c and a sharp 
decrease in the anisotropy energy. 




FIG. 2. The dependence of the anisotropic part of the piezoelectric correction to the energy of a double layer Wigner crystal 
in a (001) plane versus the direction of the Ri vector for the values of 77 specified; a, a simple hexagonal phase and a rectangular 
phase at small 77; b, a rectrangular phase at large 7? and a square phase; c, a rhombic and double hexagonal phases near the 
first-order transition point. The energy is given in units of e 2 Xv / ^/2e- 



G 



In the case of two-dimensional electron layers parallel to the (0—11) plane, the calculation gives the following results. 
A minimum of the energy of the simple hexagonal structure is reached when one of the primitive lattice vectors is 
oriented along the [100] axis. When the distance between layers is increased and the rectangular structure is formed, 
there arise two distinct equilibrium orientations, one of which corresponds to a local minimum (Ri directed along the 
[100] axis) and the other to a global minimum (Ri directed at an angle (3 ~ 60° to the [100] or [—100] axis). Near the 
points of transition to the square phase the local minimum vanishes, and a rapid reorientation of the electron lattice 
occurs. The energy of the square phase is minimum when one of the primitive lattice vectors is directed at an angle 
(3 = 45° to the [100] axis, i.e., to the vector c, which is parallel to the [100] or the [011] axis. After the transition to 
the rhombic phase the direction of the vector c parallel to the [011] axis corresponds to the global minimum, while 
the direction parallel to the [100] axis corresponds to a local minimum. As r/ is increased, the local minimum splits 
into two, corresponding to a direction of the vector c at an angle ±/3 (/3 < 30°) to this axis. At the transition to the 
double hexagonal phase all three minima become equal (the orientation of Ri changes abruptly at this transition) . 




[100] [111] [011] [-111] [-100] 



FIG. 3. The anisotropic part of the energy of the Wigner crystal in a (0-11) plane versus the direction of the Ri vector for 
the values of r/ specified; a, a simple hexagonal and a rectangular phases; b, a square and rhombic phases near the second-order 
transition point; c, a rhombic and double hexagonal phases near the first-order transition point. 
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In a bilayer structure lying parallel to the (111) plane the minimum of the energy of the simple hexagonal and 
rectangular phases correspond to a direction of Ri at an angle [3 — 30° to one of the twofold axes lying in the (111) 
plane. As the point of transition to the square phase is approached, the anisotropic contribution decreases sharply. 
A sharp reorientation occurs near the transition point. For the square phase a minimum of the energy is reached 
when the vector Ri is directed at an angle (3 = ±15° to one of the twofold axes. At the transition to the rhombic 
phase the anisotropy again increases. The energy of the rhombic phase is minimum in the case when the vector c is 
oriented along one of the twofold axes. At the transition to the double hexagonal phase the orientation of the vector 
c changes abruptly — it deviates somewhat from the twofold axis. We note that in the double hexagonal phase the 
energetically most favorable orientation of the lattice vectors differs from the case of the simple hexagonal phase. The 
latter is due to the absence of a center of inversion in the double hexagonal structure. 




FIG. 4. The anisotropic part of the energy of the Wigner crystal in a (111) plane versus the direction of the Ri vector 
for the values of r\ specified; a, a simple hexagonal and a rectangular phases; b, a rectangular and and a rhombic phases near 
the second-order transition point and a square phase; c, a rhombic and double hexagonal phases near the first-order transition 
point. 



The results show that the orientation of a bilayer Wigner crystal in a piezoelectric matrix is determined by the plane 
in which the electron layer is located, the type of electron lattice, and the value of the parameter rj. Interestingly, the 
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abrupt (jumplike) character of the reorientation is observed not only at the transition from the rhombic to the double 
hexagonal phase (this is the expected effect, since it accompanies a first-order transition). A very rapid reorientation 
also takes place at the second-order transition from the rectangular to the square phase. 

The value of the anisotropy energy is determined by the parameter %, which in GaAs is of the order of 2 x 10 -4 . The 
typical difference between the Coulomb energy in the different phases is of the order of 10~ 2 e 2 ^/n/e per electron, [f7]j|] 
i.e., according to the results reported, the piezoelectric interaction in the system under study is approximately two 
orders of magnitude smaller; hence, it has a weak influence on the phase diagram and only determines the orientation 
of the electron crystal. Nevertheless, in other systems in which the value of the piezoelectric modulus is larger, one 
can expect a radical rearrangement of the phase diagram. Analogous effects for a monolayer system were discussed 
in Ref. |jj . The approach considered in the present paper makes it possible to investigate this possibility in detail for 
the case of bilayer electron crystals. 

APPENDIX 



Let us transform to the rapidly convergent form of the expression 

:imtl>R 



s ±™ = E — 5— > ( A -l) 



R^O 

where m > 0. We introduce the function 



iq(R+r)±imi/) r +R ±imip r 

T ±m (r )q) =e-^5:^ FTSi — (A.2) 

|r+R| 1 1 



S ±m = lim T ±m (r,q). (A.3) 

r— >0,q^0 



such that 



We use the identity 

7(^,s) + r(=*M 

r(^±i) _ ' 1 ' 

With allowance for (A. 4) the quantity T± m can be written in the form a sum 

2± m (r,q) =T ±mil (r,q)+T ±m)2 (r,q), (A.5) 

where 

^ _ e *q«-±''"^r+Rr(g^,7m|r + R| 2 ) e ±m ^ j{^-,nnr 2 ) 
T ±m ,i(r,q)-2^ |r + R| p^m±i) r T(m±l) ' (A ' 6) 

jq(r+R)±imVr+R, ryfStkl nn \ r + R|2) 

W.q) = ^£ ° | f + R , (A.7) 

We note that the last term in (A. 6) vanishes in the limit r — > 0. 

For the transformation T± m 2 we substitute the definition of the function j(a, x) into (A.7): 



T ±m , 2 (r, q) = ——e-** / V |r + R \™ e in(r+n)±imi, r+R -e\r + nf 

H— ; Jo R 

Expanding (A. 8) in a Fourier series in the reciprocal lattice vectors, we get 



(A.8) 



T± m , 2 (r,q) - 7T7^T, e ~' l{Cl+G)r / / d 2 pp m exp[zp(ci + G)±zr m Pp-ep 2 }- (A.9) 

1 I - 2~ > G 



<) 



Evaluating the integral over p, we find 

|q + G| 2 



2 y G 

Making the change of variables £ = |q + G\/2t, we arrive at the form 



4£ 



2 



). (A.10) 



J 1 , ,( r a )_ r J^_y -,(q + G)rii m * 1+G 1 f ™+l |q+G| 2 

±m ' 2l ' qj " r(^)^ 6 |q+G| 1 2 ' Ann h [AAi) 

Substituting formulas (A.l), (A. 3), (A. 7), and (A. 11) into Eq. (12) and introducing the function (17), we arrive at 
Eq. (15). The term with G = in (A. 11) drops out, since it cancels with the intralayer interaction with the positive 
neutralizing background that appears in . 

An analogous transformation may be done for the sum 

^ d l \R + c\ m e ±im ^+° 

l ' ±m ~ L> (|R + CP + d2 )(i+m+ l)/2 • 

Here we also use the identity (A. 4), with m replaced by m + I. We get 

Sl,±m = 2l,±m,i(0,0)+ lim T,, ±m , 2 (r,q), (A.13) 

where 

, nn ,_^ d'lR + cre^^ r[^,^n(|R + c| 2 +d 2 )] 

Ji,±m,HU,U; -2^ (| R + c |2 +d 2)(i+m+l)/2 r( '+™ +1 ) ' { ' 

The quantity T; i ± TOi2 (r, q), after we change to summation over G and do the integration over p, reduces to the 
form 

T,, ±m , 2 (r,q) = ^_^_^ e --c- l(q+G ) r± ^ q+GdJ (iSL^l)" 

x ^^-™- 2 exp (-^d 2 - ^4^) ' (A - 15) 



Evaluation of the integral in (A. 15) gives 



N(l,m) ,. lN m +'- 2 » 



ri, ±rn ,2(r,q)-t r( i+m+1 ) 2^ e q W+G\ 4- ^0,™)+^ 



q + G|d 



where 



2 / G s=0 
;{e -| q+ G| d[r(s + 1 ) _ sgn(/ _) 7(s + 1 , /2 }] + ^i^^+G,^ + 1 ; ; 2 )}) (A 16) 



/ ± = Ja + £l±^d, (A. 17) 

2V7rn 



,m—l l — m — 2. w „, 
AT(/,m) = max(^-, ). (A.18) 

In deriving (A. 16) we took into account that in the problem considered, the parameters I and m have the same 
parity. The term with G = in (A. 16) is compensated by the interaction with the positive background of the adjacent 
layer. Using (A.13), (A. 14), and (A. 16), we can write the contribution of the interlayer interaction in the form (16). 
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